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REMARKS/ARGUMENTS 

In response to paragraph 2 of the Office Action (page 2), applicants have now removed all the new 
matters (equations, functions, matrices, etc) as the Examiner required. Applicants submit that the new 
Specification includes no new matter, and respectfully request withdrawal of the objection. 

In response to paragraph 3 of the Office Action (page 3), a marked up version of the amended 
application is attached with this submission. Applicants believe this complies with the Examiner's 
requirement. 

In response to paragraph 4 of the Office Action (page 3), applicants have added Figure 1 and Figure 2, 
so as to make the specification easier to understand. Please note that 

(1) Figure 1 and Figure 2 only show well-known concepts in automatic control; and 

(2) Although the present invention can be clearly described in words as: 

"This invention chooses the best values for the tuning parameters in a PI D controller 
or a linear controller in such a way that the largest absolute value of all poles of said 
discrete-time closed-loop transfer function from said set-point r(k) to said process 
variable v(k) i s minimized subject to. if any, user-specified constraints on one or more 
of the tuning parameters ." (see the Specification or the Claims) 

this invention cannot be described by drawings, because it contains abstract mathematical concepts 
(see the above underlined words) that cannot be described by drawings. However, anyone skilled in the 
art can easily understand the invention as described in the Specification and Claims and can easily 
carry out this invention by using well-developed optimization algorithms in mathematics or 
well-developed computer programs such as the "Optimization Toolbox" in Matlab developed by The 
MathWorks Inc. (Matlab is a registered trademark of The MathWorks Inc.). The applicants do not need 
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to teach these well-known prior arts. 

In response to paragraph 5 of the Office Action (page 4), applicants submit that the new specification is 
now in full compliance with the layout suggested by the Examiner. 

In response to paragraph 6 of the Office Action (page 8), applicants submit that the new dependant 
claims now refer to other claims in the alternative only and comply with 37 CFR 1.75 (c). Applicants 
respectfully request re-examination of the new claims and withdraw of the objections. 

In response to paragraphs 7, 8 and 9 of the Office Action (pages 8-9), applicants have removed the 
equations, transfer functions, and matrices as required by the Examiner. Anyone skilled in the art 
knows how to find the closed-loop or open-loop transfer functions, the poles of a transfer function, the 
definition of PID controllers and linear controllers and their structures, and the use of well-developed 
optimization algorithms in mathematics and optimization computer programs such as the "Optimization 
Toolbox" in Matlab developed by The MathWorks Inc. to solve the optimization problem formulated in 
the present invention (Matlab is a registered trademark of The MathWorks Inc.). The applicants do not 
need to teach these well-known prior arts . For example, if an invention needs the calculation of the 
value of sin(x), the inventor does not need to teach the definition of sin(x), neither does he need to 
teach the various mathematical methods, algorithms or computer programs for finding the value of 
sin(x), since these are all well-known prior arts. In Matlab, the function 'Tminimax" can be used to 
directly solve the optimization problem (minimax problem) formulated in this invention (see the newly 
cited book "Optimization Toolbox User's Guide" authored by Coleman et al. that comes with the 
commercial software product "Optimization Toolbox" developed by The MathWorks Inc). Matlab and the 
associated "Optimization Toolbox" are well known among people skilled in the art. Therefore, applicants 
do not need to teach how to minimize the maximum of absolute values of all poles of a discrete time 
closed-loop transfer function from a set point to a process variable. 

Therefore, applicants respectfully submit that the new claims have overcome the defects as pointed out 
by the Examiner and are in full compliance with 35 U.S.C. 112 since anyone skilled in the art can 
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implement this invention without any difficulty. Applicants respectfully request re-examination of the 
new claims and withdrawal of the rejections and criticisms as stated in paragraphs 7, 8 and 9 of the 
Office Action. 



In response to paragraph 10 of the Office Action (page 9), applicants revised the specification so that it 
now relies on the prior art of record cited. 



Applicants have read through the following new reference cited by the Examiner: 

A. Cichocki et al., "Neural Networks for Solving Systems of Linear Equations II: Minimax and 
Least Absolute Value Problems", IEEE Transactions on Circuits and Systems II: Analog and 
Digital Signal Processing, Volume 39, Issue 9, Sept. 1992, Pages 619-633. 

A careful study of this prior art shows it deals with the problem of solving a system of linear equations 
by minimizing the largest absolute error using the neural networks technology. It is not apparent that 
this method can be used to solve the minimax problem formulated in this invention since the minimax 
problem in this invention is totally different. However, there are many well-developed optimization 
methods in mathematics that can directly solve the minimax problem formulated in this invention (see, 
e.g., the many cited references mentioned in the "DETAILED DESCRIPTION OF THE INVENTION" 
section). Indeed many computer software packages such as the "Optimization Toolbox" in Matlab 
developed by The MathWorks Inc. can solve the minimax problem formulated in this invention without 
any difficulty. Matlab and its various toolboxes such as the "Optimization Toolbox" are very popular 
among people skilled in the art. All these well-developed optimization methods are well-known prior art 
and do not need to be taught in this invention. 



Applicants respectfully submit that the new specification now relies on the cited references and 
therefore complies with the requirements in paragraph 10 of the Office Action. 
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Conclusion 



In order to comply with the Office Actions, Applicants have re-written the specification and claims. For 
all the aforesaid reasons applicants submit that the new specification and claims do not contain new 
matters. Drawings are added. The layout of the specification is now as suggested by the Examiner. The 
specification and claims enable anyone skilled in the art to carry out the invention without any difficulty. 
Applicants respectfully submit that this application is now in full condition for allowance. 

However, applicants are willing to accept any suggestions that the Examiner believes necessary in 
order to avoid further proceedings. 

Very respectfully, 





Danyang Liu 



Wuyang Liu 



#108, 10615-111 Street 



P.O. Box 112 



Edmonton, Alberta T5H 3G1 



316, San Jing Lu 



Canada 



Nanchang, P. R. China 



Telephone: (780) 424-5074 
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VERSION WITH MARKINGS TO SHOW CHANGES MADE 
In the Specification: 

Please cancel the old specification in record and replace it with the new Specification (underlined) as 
follows: 

Doscr i pt i on 

T I TLE OF THE I NVENTION 

Mothods of Dosigning Opt i ma l PID Controllorc 

TECHN I CAL F I ELD OF THE I NVENT I ON 

This i nvont i on rolatos to tho dosign of the structure of a mu l tivar i ab l o P I D controllor and tho opt i ma l 
cho i co of its P I D paramotors. 

BACKGROUND OF THE I NVENT I ON 

A trad i tiona l P I D control l or is ucod to control an i ndustria l procoss. Tho procoss variab l o (PV) goos i nto 
tho PLD controllor, wh i ch ca l culator tho contro ll or output (CO) accord i ng to a P I D control oquation. Th i s 
CO is th o n convert e d to an a na l og sign al , wh i ch i s sent to the process so that tho s ai d PV can track a 
usor opoc i fiod va l uo callod cot po i nt (SP). Tho sa i d SP can chango w i th timo. Tho porformanco of a P I D 
contro l lor dopond c on tho cho i co of i ts thr o o PID paramotors. For indopondont form of PID contro l lors 
the s e throo P I D param e t e r s aro tho proport i ona l ga i n Kp, th e int e gra l ga i n K i , and tho d e r i vat i v e gain 
Kd. For d e p e ndent form of PID contro lle rs thoso thro o P I D p a ramet e rs a ro tho ga i n K, int e gral t i m o T i , 
and dorivativo t i mo Td. I n trad i t i ona l PID contro l lors tho sa i d PV, SP, CO, and PEP paramotors aro a ll 
sca l ars. Wo call th i s k i nd of P I D controllors tho s i nglo - input singlo - output (SISQ) P I D contro ll or c . Tho 
Zi e g le r - N i chols P I D contro l l o r tun i ng method is tho major on o of tho many m e thods for find i ng th o 
va l ues of P I D paramotors. 

DETA I LED DESCR I PTION OF THE INVENT I ON 
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I n th i s i nvontion tho S I SO P I P contro ll or i s oxtondod to tho multiplo -i nput mu l t i plo - output (M l MO) P I P 
contro ll or that has n procooo var i ablos PV1, PV2, . . . , and PVn and m contro l lor outputc CQ1, 
CQ2, . . . , and COm t whoro m and n aro pos i t i vo i ntogors, Correspond i ng to PV1 , PV2, .... and PVn 

th e r e ar o n s o t points SP1 , SP2 and SPn. In th i s cas e PV b e com e s a v e ctor w i th PV1, PV2 , 

and PVn being i ts first, second, . . . , and n - th component, CO becomes a voctor with CQ1 , CQ2, 
and CQm boing i ts first, socond, . . . , and m-th Gomponont, SP bocomos a voctor with SP I , SP2, 
and SPn being its f i rst, socond, . . . , and n - th componont, and tho P I P contro l oquat i on bocomoc 
CO(k)-CO(k 1)^KrSP(k)*T+Kr a (k ) 1)+K2* a (k,2) +. . . +Kj*a(k,j), wh o ro k is th o discrete tim e , T i s tho 
s a mp li ng per i od, j i s a posit i vo i ntogor, K1, K2, . . . , Kj are m by n P I P paramotors, a(k,1)-[ PV(k)]*T, 
and a(k,j)-[a(k,j - 1) a(k 1 j - 1))/Tfor f> or-2. I t is important to noto that 

1. an M I MO P I P contro l lor is a b l e to tak e into account th o interaction a mong the n process v a riab lo s 
and m contro ll or outputs, which can not bo achiovod by simp l y apply i ng S I SO P I P control l ers to oach of 
tho n procoss variab l os, and 

2. there i s no s o t po i nt i n any of a(k,1), a(k,2), .... and a(kj), wh i ch can avoid th o unwant e d sudden 
change i n CO when SP changes with t i me. 

Tho noxt prob l om of dos i gning tho optima l PIP controllor i s to find tho bost valuos for tho P I P 
paramot o rs K1 , K2, . . . , and Kj . An optim i zat i on based method for so l v i ng this prob l em cons i sts of tho 
follow i ng four stops: 

1 . Convort tho P I P contro l oquat i on i nto d i scroto tim e form i f i t i s not i n discr o to timo form. 

2. Bu i ld a d i screte t i m e li n e ar mode l for th o proc o ss th a t i s to b o controll e d by th o s a id PIP control le r. - 

3. Form tho discroto t i mo c l osod loop transfor function from sa i d voctor SP to said voGtor PV. - 

4. F i nd tho bost P I P paramotors by us i ng an opt i m i zat i on a l gor i thm wh i ch m i n i mizos tho l argoct 
modulus of a l l polos of th e d i scr e t e t i m e c l osed loop transf o r funct i on obta i n e d at stop 3, whoro tho 
modu l us of a po l o i s defined to bo tho abso l ute va l ue of tho complex number wh i ch represents tho polo. 
If tho P I P param e ters aro subjoct to somo constra i nts, thon a constrainod opt i mizat i on a l gorithm can bo 
usod which m i nim i zos tho largost modu l us of all po l o s of tho d i scroto timo closod l oop transfor funct i on 
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obta i nod at stop 3 and at tho samo t i mo guarantees that a l l u c or proscr i bod constra i nt s on tho PID 
paramotorc aro caticfiod. 

P I D contro lle rs with the i r param e t e rs so obtain o d guarante e that PV can track SP qu i ck l y 

SPECIFICATION 

TITLE OF THE INVENTION 

Methods of Designing Optimal Linear Controllers 

CROSS-REFERENCE TO RELATED APPLICATIONS 

Domestic priority data: This application is a 371 of PCT/IB01/01002 06/07/2001 
Foreign applications: June 20. 2000 fCAl Canada 2.311 ,268 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT 
Not Applicable. 

INCORPORATION-BY-REFERENCE OF MATERIAL SUBMITTED ON A COMPACT DISC 
Not Applicable 

BACKGROUND OF THE INVENTION 
Technical Field of the Invention 

This invention relates to the design of optimal linear controllers and PID controllers. 
Description of Related Art 

Figure 1 shows a process 1 controlled by a PID (proportional-integral-derivative) or a linear 
controller 2. The PID controller or linear controller 2 means a linear system, usually implemented in a 
computer or as an electronic circuit, that receives the process variable signal v(k) and set-point signal 
r(k) as its inputs and calculates the controller output signal u(k) according to a PID control eguation or a 
linear control eguation, where k is the integer discrete time variable. For a multiple-input and 
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multiple-output (MIMO) process, the process variable y(k) is an n-dimensional vector with each of its 
components being a scalar process variable and the controller output u(k) is an m-dimensional vector 
variable with each of its components being a scalar controller output, where n and m are positive 
integers. It is desired that the performance of the controller should be such that, after the controller 
output signal u(k) is sent to the process 1. the process variable y(k) should approach the set-point rlk) 
quickly and smoothly as the time variable k increases, where r(k) is an n-dimensional set-point signal 
(also known as the reference signal or the command signal, etc.) 

There are many types of PIP controllers, depending on the use of different types of PIP control 
equations. All types of PIP controllers can be viewed as special cases of linear controllers. By definition, 
a discrete-time linear controller means a controller whose output u(k) is a linear function of past 
controller outputs, current and past process variables, and current and past set-point variables, as 
shown in Figure 2. This definition is well known to any person skilled in the art (see, e.g.. Figure 2.4 and 
Assumption 2.3 on page 29 of the cited book "Linear Controller Design" authored by Boyd and Barratt 
and published by Prentice Hall in 1991). 

Once the structure of its control eguation is properly selected, the performance of a PIP controller 
or a linear controller depends mainly on the choice of coefficients in the control eguation. The 
adjustable numbers in the coefficients are called tuning parameters (also known as tuning coefficients, 
tuning gains or tuning constants, etc.) 

How to properly choose tuning parameters for a PIP controller is a problem that has attracted a lot 
of studies ever since PIP controllers became widely used in industry in the early 1940s. The 
Ziegler-Nichols tuning methods developed by Ziegler and Nichols in 1942 (see the Ziegler and Nichols 
reference), either in their original form or in some modification, are still widely used in industry. Other 
model-based optimization methods choose the tuning parameters by minimizing some well-known 
control performance index such as the integrated absolute errors (IAE), the integrated sguared errors 
(ISE). the integrated time absolute error (ITAE). etc. (see US patent 5.453,925). However, practice 
shows that all these methods often lead to the undesired oscillatory control results. 

Many tuning methods are based on the continuous-time transfer function analysis (see, for 
example, the cited US patents 6.434.436 B1 to Adamv et al.. 5.866,861 to Raiamani et al.. 4.539.633 to 
Shigemasa et al.. 4.563.734 to Mori et al.. the cited research papers by Celentano et aL Zhuanq et al.. 



PACE 14/33 * RCVD AT 7/23/2004 11:1 0:57 AM [Eastern Daylight TimeJ* SVR:USPTO-EFXRF-1/0" DNlS:8729306 * CSID:7804245074 * DURATION (mm-ss): 17-00 



07/23/2004 09:10 



Dan Liu 



PAGE 15/33 



Application/Control Number 10/069,364, Art Unit 2121 Page 5 

Dorf et al.. Katebi et aL and Saeki et al.. and the cited book "Linear Controller Design" authored bv 
Bovd et al.). Therefore these methods cannot guarantee the optimal performance of controllers in the 
discrete-time domain. This is a real problem since nowadays almost all PIP or linear controllers are 
implemented in digital computers and therefore they are working in the discrete-time domain, not the 
continuous-time domain. 

Other tuning methods in discrete-time domain cannot guarantee that the largest absolute value of 
all the poles of the closed-loop transfer function is minimized. Therefore they cannot guarantee the 
optimal performance of the controller (see, e.g.. the cited US patent No. 5,680,304 to Wang et al. and 
the cited research paper bv Yamamoto et aL. etc.). 

BRIEF SUMMARY OF THE INVENTION 

This invention chooses the tuning parameters in a PIP controller or a linear controller by minimizing 
the maximum of absolute values of all poles of the discrete-time closed-loop transfer function from said 
set-point r(k) to said process variable v(k) subject to, if any, user-specified constraints on one or more 
of the tuning parameters. When the tuning parameters are chosen this way, the PIP controller or linear 
controller can guarantee that the process variable v(k) tracks the set-point r(k) smoothly and ouicklv as 
time k increases. 

BRIEF PESCRIPTION OF THE PRAWINGS 

Figure 1 shows a process 1 controlled by a PIP controller or a linear controller 2. The PIP controller 
or linear controller is in discrete-time form. 

Figure 2 shows a PIP or a linear controller. Bv definition, the controller output u(k) of a linear 

controller is a linear function of u(k-1), u(k-2) u(k-a). y(k). v(k-1) y(k-b), r(k). r(k-1^ r(k-2) 

r(k-c). wherein a>0, b and c are integers (This definition is well known to any person skilled in the art 
and is given in Figure 2.4 and Assumption 2.3 on page 29 of the cited book "Linear Controller Pesian" 
authored bv Boyd et al.) 

PETAILED PESCRIPTION OF THE INVENTION 
From now on it is always assumed that: 
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(1) The open-loop discrete-time transfer function of the process 1 is known, and 

(2) The structure of the control equation of the PIP controller or the linear controller has already 
been properly chosen and is already in discrete-time form. 

It is then easy for anyone skilled in the art to find the discrete-time closed-loop transfer function from 
the set-point r(k) to the process variable v(k). as shown in Figure 1. This invention chooses the best 
values for the tuning parameters in a PIP controller or a linear controller in such a way that the largest 
absolute value of all poles of said discrete-time closed-loop transfer function from said set-point r(k) to 
said process variable v(k) is minimized subject to. if any, user-specified constraints on one or more of 
the tuning parameters. This choice guarantees that the process variable v(k) tracks the set-point rfk) 
smoothly and guickly as time k increases. 

The above description contains a minimax optimization problem. The general minimax problem 
has been very well studied and many successful numerical methods and algorithms have been 
developed, see, e.g., the cited publications by Charalambous et al.. Conn et al., Conn. Pillo, Gigola et 
aL Polvak. Polak. Polak et al.. Vardi, Zang. Murray et al.. Kaufman et aL Reemsten. Zhou et aL and 
Laskari et al. Successful commercial computer programs such as the "Optimization Toolbox for use 
with Matlab" developed by The MathWorks Inc. can directly be used to solve the minimax problem as 
formulated in this invention without any difficulty (see the cited book "Optimization Toolbox User's 
Guide" authored by Coleman et al. and published by The MathWorks Inc.). The "Optimization Toolbox" 
and Matlab have been well known among people skilled in the art. It is easy for anyone skilled in the art 
to solve the minimax problem directly using the "Optimization Toolbox", the methods in the cited 
publications mentioned above, or any other method. 

In the Claims: 

Please cancel all claims in record and replace them with the new claims 22-33 (underlined) as follows: 
Wq c l a i m: 

1. AM I MO (mu l t i p l o - input mu l t i p l o - output) P I P contro ll or which has an n - d i monsiona l procoss variab l o 
voctor PV with tho n process var i ab l es FV1, PV2, .... and PVn boing i ts first, socond, . . . , and n th 
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componont, an n - d i mons i onal sot po i nt voctor SP w i th tho n s ot po i nts SP1 f SP2, . . . , and SPn bo i ng 
its f i rst, oocond, . . . , and n th componont, and an m - d i monsiona l contro l lor output voctor CO with tho m 
contro ll or outputs CQ1 , CQ2, . . . , and COm bo i ng it s first, socond, . . . , and m-th componont, whom m 

and n are pos i t i v e intogors, and m wh i ch — ±be P\D contro l e quation is 

CO(k)-CO(k I^KrSPQcrT^Kra^lH^aCk^)* - . . , +Kj*a(kj), whoro k is tho discrete t i mo, T i s tho 
samp li ng por i od, j i s a pos i tivo i ntogor, K1, K2, . . . , Kj aro m by n P I D paramotors, a(k,1)-[ - PV(k)]*T, 
and a(kj)-[a(k,j - 1) - a(k - 1 j - 1)]/Tfor or "2. 

2. An M I MQ P I D contro ll or of cla i m 1, in which tho m by n PID paramotors K1, K2, . . . , and Kj aro 
obta i nod by us i ng an optimizat i on a l gor i thm wh i ch m i n i mizos tho largo s t modu l us of all po l os of tho 
di s cr e t e t i m e c l o s ed l oop transf e r funct i on from said SP to s a i d PV. 

3. An MIMQ P I D contro ll or of c l a i m 2, whoro i n tho said opt i m i zat i on algor i thm i s a conctrainod 
opt i m i zation algorithm wh i ch m i n i m i zos tho l argost modu l us of al l po l os of tho d i scroto timo c l osod loop 
transf e r funct i on from sa i d SP to sa i d PV and a t th e sam e ti me guarant ee s that th e usor proscribed 
constraints on tho PID parameters aro sat i sfied. 

4. An M I MQ P I D contro ll or of cla i m 1, whoroin somo or all of tho torms K2*a(k,2), K3*a(k,3), . . . , and 
Kj* a (k,j) that app e ar on tho r i ght - hand s i d e of th e PID contro l e quat i on ar e r e mov e d, for e xamp l e, a P I D 
control l or with — its contro l oquat i on — b ei ng CQ(k)=CQ(k - 1)^K1*SP(k)*T^Kra( — k,1)-CO(k 1) 
*K1*[SP(k) - PV(k)]*T, wh i ch i s als o ca ll od a 1 - on l y control l or, and a P I D contro ll or with its contro l 
oquat i on — being — CO(k)-CO(k - 1)+KrSP(k)*T+Kra(k,1)+K2*a(k 1 2)-CO(k - 1)+K1*[SP(k) PV(k)]*T K 
2*[PV(k) PV(k - 1)], wh i ch i s al s o c allod a P I control lo r, e tc . 

5. An M I MQ P I D controllor of c l aim 4, whoro i n tho rema i n i ng PID paramotors aro obtain o d by us i ng an 
optim i zation algor i thm wh i ch m i nim i z e s th e l arg e st modu l us of a ll pol e s of the d i scroto timo clos o d loop 
transfer funct i on from sa i d SP to sa i d PV. 

6. An M I MQ P I D controHor of c l a i m 5, whoroin tho sa i d opt i m i zation a l gor i thm is a constrainod 
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optim i zat i on a l gorithm wh i ch m i nimizos tho l argo s t modulus of al l po l os of tho d i scroto timo c l osod l oop 
transfor function from said SP to said PV and at tho camo t i mo guarantoos that tho ucor proccriboc t 
con s tra i nts on tho P I D p a rameters are sat i sfiod. 

7. A P I D control l er of cla i m 1, whoroin said PV, said SP, sa i d CO, and sa i d PID parameters aro al l 
s ca l ars, a nd m-n-1 . 

8. A P I D contro ll er of c l a i m 2, whoro i n sa i d PV, s aid SP, sa i d CO, and sa i d P I D paramet e rs a r o a ll 
sca l ar s , and m-n~1. 

9. A P I D contro ll or of c l aim 3, whoro i n said PV, sa i d SP, said CO, and s a i d P I D p a ram e t e r s aro a ll 
scalars, a nd m=n=1. 

10. A P I D contro ll or of c l aim A t whoro i n sa i d PV, sa i d SP, said CO, and sai d P I D paramotors aro a ll 
sca l ars, and m=n = 1. 

11. A P I D contro ll or of cla i m 5, whoroin said PV, said SP, said CO, and sa i d P I D paramotors aro al l 
s calars, and m~n-1 . 

12. A P I D contro l lor of cla i m 6, whoro i n sa i d PV, said SP, caid CO, and sa i d PID paramotors aro a l l 
s c al ar s , and m-n-1. 

13. A method of find i ng tho optimal P I D parameters or tun i ng constants for a ny P I D contro l l o r or a li n e ar 
contro l lor by m i n i miz i ng tho l argost modu l u s of a l l po l os of tho d i scroto - t i mo c l osod -l oop transfor 
funct i on from so t po i nt SP to process variab l o PV. 

14. A method of C l aim 13, whoroin tho opt i ma l P I D p a r a m e t e rs or tun i ng constants aro obtained by 
m i nimiz i ng tho l argost modu l us of tho discroto - t i mo closod -l oop transfor funct i on from said SP to ca i d 
PV subjoct to tho con c tra i nt that somo or a ll of tho P I D paramotors or tuning constant s aro within ucor 
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s poc i fiod admissible * rangos. 

15. A P I D contro l lor of C l aim 1 or a ny PID contro l lor or l i noar contro l lor, whoro i n tho P I D paramotor s or 
tuning constants aro obta i ned such that 

(1) tho maximum roal part of a ll po l os of tho continuous - t i me c l osod -l oop transfer funct i on from sa i d SP 
to said PV is m i n i m i zod, pos si b l y s ub jo ct to any usor spocifiod constraints on tho PID paramotors or 
tuning constants, s uch as somo or a ll of olomonts i n tho PID paramotors or tun i ng constants shou l d 
bo w i th i n usor specified rangos, whor o th o s ai d P ID par a m ete rs or tuning constants can b o matr i ces 
or sca l ars, or 

(2) tho max i mum magn i tud o of s omo or a ll olomonts i n tho PID paramotors or tun i ng constants i s 
m i n i m i z e d, subjoct to the con s tra i nt that th e maximum r eal part of a ll po l os of tho sai d 
cont i nuous - t i m e c l os e d - loop transf e r funct i on from said SP to said PV is not la rg e r than a usor 
spoc i fiod numbor. 



We claim: 

22. A method for determining the optimal tuning parameters in a linear controller, wherein 

1) said controller receives an n-dimensional process variable signal v(k) from a process and an 
n-dimensional set-point signal r(k), calculates an m-dimensionat controller output u(k) 
according to a linear control equation, and sends said u(k) to said process, where k is the 
integer discrete time variable and n and m are positive integers: 

2) said tuning parameters are the adjustable numbers in the coefficients in said linear control 
eguation that are to be determined: and 

3) said method finds the optimal values for said tuning parameters by minimizing the maximum of 
absolute values of all poles of the discrete-time closed-loop transfer function from said 
set-point r(k) to said process variable vfk). 

23. A method as in Claim 22. wherein said minimization of the maximum of absolute values of all poles 
of said discrete-time closed-loop transfer function is subject to user-specified constraints placed on 
one or more of said tuning parameters. 

24. A method as in Claim 22. wherein said controller output ufk) = u(k-1) ± K i *r(k)*T ± K i *a(k,1) + 
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K?*a(k.2) + ... + Kp*a(k,p), wherein k is the discrete time variable. * is the multiplication operator. T 

is the sampling period, p is a positive integer, the m by n matrices K i , K? and K p are tuning 

parameters. a(kA) = r-yfl01*T. and afk. p) = fafk.p-D - a(k-1,p-1)1 / Tfor p> or =2. 

25. A method as in Claim 23, wherein said controller output u(k) = u(k-1) + K i *r(k)*T ± K i *a(k,1) + 
K ? *a(k,2) + ... +K P *a(k.p). wherein k is the discrete time variable. * is the multiplication operator. T 

is the sampling period, p is a positive integer the m by n matrices K i . K? and K P are tuning 

parameters. a(k.1) = f-v(k)rT. and a(k. p) = fa(k,p-1) - a(k-1 ,p-1)1 / T for p> or =2. 

26. A method as in Claim 22. wherein said linear controller is a PIP (proportional-integral- derivative) 
controller. 

27. A method as in Claim 23. wherein said linear controller is a PIP controller. 

28. A linear controller as in Claim 22 with its tuning parameters determined therein. 

29. A linear controller as in Claim 23 with its tuning parameters determined therein. 

30. A linear controller as in Claim 24 with its tuning parameters determined therein. 

31. A linear controller as in Claim 25 with its tuning parameters determined therein. 

32. A PIP controller as in Claim 26 with its tuning parameters determined therein. 

33. A PIP controller as in Claim 27 with its tuning parameters determined therein. 

In the Abstract 

Please cancel the old Abstract in record and replace it with the new Abstract (underlined) as follows: 

Mothod s of do si gn i ng tho structuro of mu l t i p l o -i nput mult i p l o-output (MIMO) P I P contro ll ers and 
m e thods of find i ng th e optimal va l u es for tho M I MO P I P param o t e rs aro d is c l osod. Th e optima l va l uoc 
of M I MO P I P par a m e t e rs a r e obt a in e d b y using an opt i mizat i on a lgor i thm which m i n i miz es th o larg e st 
modu l us of al l po l os or tho discroto t i mo c l osod loop transfor function from sot point SP to procooc 
variab l o PV, with or w i thout usor proscr i bod con s tra i nts on tho P I P paramotors. Mothods of dosign i ng 
th e structur o of s i ng le i nput sing l o - output (S I SQ) P I P control l ers and m o thods of find i ng t i l o opt i ma l 
va l ues for S I SQ PIP paramotors aro a l so disclosod as spec i a l M I MO P I P contro lle r cases. 

Methods of designing optimal discrete-time PIP (proportional-integral-derivative) controllers and linear 
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controllers are disclosed. The optimal values of the tuning parameters in a PIP controller or a linear 
controller are determined by minimizing the maximum of absolute values of all poles of the 
discrete-time closed-loop transfer function from the set-point to the process variable subject to, if any, 
user-specified constraints on one or more of the tuning parameters. 

In the drawings: 

Please add the attached Figure 1 and Figure 2. 
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